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ABSTRACT 

The early evolution of Earth's atmosphere and the origin of life took place at a time when physical 
conditions at the Earth where radically different from its present state. The radiative input from the 
Sun was much enhanced in the high-energy spectral domain, and in order to model early planetary 
atmospheres in detail, a knowledge of the solar radiative input is needed. We present an investigation 
of the atmospheric parameters, state of evolution and high-energy fluxes of the nearby star k 1 Cet, 
previously thought to have properties resembling those of the early Sun. Atmospheric parameters 
were derived from the excitation/ionization equilibrium of Fe I and Fe II, profile fitting of Ha and 
the spectral energy distribution. The UV irradiance was derived from FUSE and HST data, and the 
absolute chromospheric flux from the Ha line core. From careful spectral analysis and the comparison 
of different methods we propose for k 1 Cet the following atmospheric parameters: T e g = 5665 ± 
30 K (Ho profile and energy distribution), \ogg — 4.49 ± 0.05 dex (evolutionary and spectroscopic) 
and [Fe/H] = +0.10 ± 0.05 (Fe II lines). The UV radiative properties of k 1 Cet indicate that its 
flux is some 35% lower than the current Sun's between 210 and 300 nm, it matches the Sun's at 
170 nm and increases to at least 2-7 times higher than the Sun's between 110 and 140 nm. The 
use of several indicators ascribes an age to k 1 Cet in the interval ~ 0.4-0.8 Gyr and the analysis 
of the theoretical HR diagram suggests a mass ~ 1.04 M . This star is thus a very close analog 
of the Sun when life arose on Earth and Mars is thought to have lost its surface bodies of liquid 
water. Photochemical models indicate that the enhanced UV emission leads to a significant increase 
in photodissociation rates compared with those commonly assumed of the early Earth. Our results 
show that reliable calculations of the chemical composition of early planetary atmospheres need to 
account for the stronger solar photodissociating UV irradiation. 

Subject headings: stars: chromospheric activity — stars: abundances - stars: late-type - techniques: 
spectroscopic — planets: atmospheres 



1. INTRODUCTION 

The irradiation from the parent star is, by far, the most 
important source of energy in planetary atmospheres. 
Most of the physical and chemical properties of the at- 
mosphere of a planet are largely driven by the stellar in- 
put, which ultimately determines its structure, composi- 
tion and, even, its mere existence. Most of the radiation 
emitted by the Sun comes from the photosphere and, 
at the solar effective temperature today (T c gQ = 5780 
K), it is the dominant source at wavelengths above 
170 nm. The solar photospheric flux is quite stable 
over short timescales (decades) and it only suffers varia- 
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tions driven by sunspots and fa culae with an amplitud e 
below 0.2-0.3% peak to peak (|Frohlich fe Lean! I2004f ). 
Over long timescales, the variations have been larger 
and are related to the nuclear evolution of the Sun. 
Model predictions indicate that the young ZAMS Sun 
could have had a luminosity some 35% lower than to- 
day (unless a scenario of heav y mass loss is considered; 
ISackmann fe Boothrovdll2003t ) . 

The energetic end of the solar spectrum (i.e., be- 
low 170 nm) is dominated by the emissions from high- 
temperature plasma in the chromosphere, transition re- 
gion and corona. Such high-energy fluxes are strongly 
variable over short- and mid-term timescales because of 
flare events, rotationa l modulation, activity cycles, etc. 
(e.g., iLean et all 119971). Over longer timescales, several 
past studies (|Zahnle fe WalkerHl982t IAvreslll997ft estab- 
lished that the young Sun's high-energy emissions were 
possibly up to several orders of magnitude stronger than 
currently. 

The Sun in Tim e programme (jDorren fc Guinanl[r99l 
iRibas et"aT1l2005L hereafter Paper I) focused on a small 
sample of carefully-selected and well-studied stellar prox- 
ies that represent key stages in the evolution of the Sun. 
This approach allowed the study in the X-ray, EUV, and 
FUV domains, where the variations are of one order of 
magnitude or more. However, the stellar proxy tech- 
nique becomes increasingly uncertain in the UV since the 
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stars are not perfect matches to the Sun (their masses are 
within 10% of 1 M©) and the expected flux variations are 
in the range of tens of a percent. For this reason, most of 
the Sun in Time stars employed in Paper I cannot be reli- 
ably used to infer the solar flux evolution at wavelengths 
longer than about 130 nm. But at the same time, the 
wavelength interval between 130 and 200 nm is an essen- 
tial energy input to planetary atmospheres since it drives 
most of the photochemical reactions (e.g.. iCanuto et ail 
119821 119831 ) . In addition, its impact on living organisms 
must also be considered, although it is unlikely that pho- 
tons in this r ange could penetrate a dense atmosphere 
(|Cockellll2000T) . 

Among the solar proxies studied in the Sun in Time, 
k 1 Cet (HD 20630, HIP 15457, d = 9.16 pc, V = 4.84) 
stands out as potentially having a mass very close to 
solar and a young age (see Paper I). This could be a 
very good analog of the Sun at the critical time when 
life is thought to have originated o n Earth 3.8 Gyr ago , 
at the start of the Archean epoch (jMoizsis et alJ ll996). 
It is also about the time when Mars lost its liquid wa- 
ter inventory at the end of the Noachian epoch some 
3.7 Gyr ago (jJakoskv fe Phillips! 120011) . For these rea- 
sons, k 1 Cet deserves attention as a possible precise 
match to the young Sun, thus providing information on 
the radiation environment that determined the proper- 
ties and chemical compositio n of the planetary atmo- 
spheres. iCnossen et al.l (|2007l) presented a study based 
on k 1 Cet with a goal set on assessing the biological im- 
plications of the high-energy radiations. The authors did 
not use a real high-energy spectrum but generated syn- 
thetic data using plasma models from an inferred emis- 
sion measure distribution. 

In this paper we carry out an in-depth analysis of 
K Cet, including its radiative properties, chemical abun- 
dances, atmospheric parameters and state of evolution, 
with the ultimate goal of understanding the Sun's UV 
emissions at a critical time in the past. The paper is 
organized as follows. In §2 we describe the spectroscopic 
analysis and derivation of atmospheric parameters, along 
with its Li abundance. The high-energy radiative pro- 
cesses and magnetic activity of k 1 Cet are discussed in 
the context of active stars in §3. In §4 we determine 
its mass and state of evolution, and in §5 we employ its 
properties to discuss the young Sun in the frame of pho- 
todissociation calculations. Our conclusions are drawn 
in §6. 

2. SPECTROSCOPIC ATMOSPHERIC PARAMETERS 
2.1. Visible spectra 

Visible spectroscopic observations were performed with 
two setups. The Coude spectrograph of the 1.60- 
m telescope of Observatorio do Pico dos Dias (OPD, 
Brasopolis, Brazil), operated by Laboratorio Nacional de 
Astrofisica (LNA/CNPq), was used to obtain spectra of 
k 1 Cet and the Sun (represented by Moon spectra) in a 
150 A spectral range, centered in the Ha line, with R 
= 20 000 and S/N ratios in excess of 300. k 1 Cet was 
observed in 2008, and the Moon spectra were acquired 
in a series of runs between 1994 and 2004, to ensure the 
absence of long-term systematics. 

Additionally, the FEROS spectrograph (Kauf er et al.l 
1999) was used in 2000 to obtain spectra of k 1 Cet and 



Ganymede (as a solar template), at a nominal resolution 
of R = 45 000. The S/N ratio per resolution element 
was 800 and 1000, respectively, for k 1 Cet and the Sun, 
between AA5600-6960, being about half as much for the 
AA4500-5500 range. As n 1 Cet is very solar-like in its 
properties, the Sun is the natural choice as the standard 
star of a differential analysis. In this approach, system- 
atic errors in line measurement, atmospheric modelling 
and the possible presence of Non-Local Thermodynamic 
Equilibrium (NLTE) effects tend to be minimized. 

Data reduction was carried out by standard procedures 
using IRAFfl Extreme care was taken in the continuum 
normalization to guarantee equivalent width (hereafter 
W\s ) measurements as free from systematic effects as 
possible. This meant chosing specific wavelength inter- 
vals by comparing, window by window, k 1 Cet with the 
solar template (Ga nymede) and the Solar Flux Atlas 
(|Kurucz et al.lll984D using t he Utrecht spectral line com- 
pilation (jMoore et al.l 1196(f ). For each window, the line 
strengths were checked to see if the wings of neighboring 
lines might be causing different depressions in the spec- 
tra. All windows considered suspect were not used. The 
remaining windows were then adjusted with low order 
polynomials. The spectral windows were always simul- 
taneously normalized (Sun and n 1 Cet) to minimize the 
possibility of errors. The W\s of Fe I and Fe II lines were 
measured by Gaussian fits. 

A series of sanity checks were performed on the mea- 
sured Was. Saturated lines were eliminated by a 2a clip- 
ping on the relation of reduced width W\/X with line 
depth, and no lines were measured beyond the linearity 
limit. Also, no trend is expected in the relation of the 
line full-width-half-maximum and reduced width, since 
the line widths are essentially set by the instrumental 
profile. For stronger lines the inability of gaussian fits 
to correctly reproduce the line wings leads to an arti- 
ficial increase of the line FWHM with reduced width, 
and no lines were measured beyond this limit. The fi- 
nal line selection comprised 90 Fe I and 12 Fe II lines. 
We follow the same proce dure discussed in detail by 
iPorto de Mello et al.l (|2008[) . The measured W A s were 
corrected by +3.6% to bring them onto a system com- 
patibl e with the Voigt- fitted solar Was of IMevlan et al.1 
(|1993f ). The reason for the correction is both a low level 
of scattered light in the spectrograph (< 2%) and also 
the inability of Gaussian fits to fully represent observed 
line profiles. 

2.2. Spectroscopic analysis 

We derived solar gf-values for the Fe I and Fe II spec- 
tral lines from a LTE, 1-D, homogeneous and plane- 
parallel solar model atmosphere from the NMAR CS grid 
i|Gustafsson et al.l I2008t lEdvardsson et al.l I1993D . The 
adopted parameters for the Sun were T c g — 5780 K, 
\ogg = 4.44, [Fe/H] = +0.00 and & = 1.00 km s" 1 , 
and we employed the Was measured off t he Ganymede 
spectr a, corrected to the Voigt scale of IMevlan et al.l 
(1 19931). The adop t ed sol ar absolute abundances are those 
of lAsplund eTafl (p009h . 

7 Image Reduction and Analysis Facility (IRAF) is distributed 
by the National Optical Astronomical Observatories (NOAO), 
which is operated by the Association of Universities for Research in 
Astronomy (AURA), Inc., under contract to the National Science 
Foundation (NSF). 
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TABLE 1 
Atmospheric parameters of k 1 Cet 



Authors 



logg [Fe/H] 



Mcthod a 



Cayrel de Strobel 
& Bentolila (1989) 
Pasqmni_et i ^il ii (1994) 
Ottman et al. (19981 
Gaidos fc Gonzalez (2002) 
Barklem et al. (2002) 
Heiter & Luck (20031 
Allende-Prieto et al. (2004) 
Luck fc Heiter (2006) 
This work 
This work 
This work 

ADOPTED 



5630 


4.50 


+0.04 


phot., Ha 


5675 


4.35 


-0.01 


phot. 


5680 






Ha, H/3 


5747 


4.53 


+0.11 


exc. & ion. 


5710 






Ha,H/3 


5750 


4.55 


+0.05 


line depths 


5564 


4.52 


-0.06 


phot. 


5700 


4.55 


+0.05 


exc. & ion. 


5780 


4.48 


+0.07 


exc. & ion. 


5685 






phot. 


5645 






Ha 


5665 


4.49 


+0.10 


sec text 



a Method: "phot." stands for photometric derivations of T G g , "Ha" 
or "H/3"; for the fitting of theoretical profiles to the Balmer lines, 
"line depths" for the calibration of T c g and line depth ratios of metal 
lines, and "exc. & ion." for the satisfaction of the excitation and 
ionization equilibrium with Fe I and Fe II lines. 



The spectroscopic atmospheric parameters of 
k 1 Cet were determined by the simultaneous satis- 
faction of the excitation and ionization equilibria of 
Fe I and Fe II. T c g was obtained by forcing the Fe I 
line abundances to be independent of their excitation 
potential. Surface gravity was determined by forcing the 
lines of Fe I and Fe II to agree to the same abundance. 
The microturbulcnce velocity was set by forcing the 
lines of Fe I to be independent of their W\s. The Fe 
abundance [Fe/H] (we use throughout the notation 
[A/B] = \ogN(A)/N(B) stai -\ogN(A)/N(B)sun, where 
N denotes the number abundance) is automatically 
obtained at the end of the iteration, and this solution 
is unique for a given set of gf values, W\s and model 
atmospheres. Here we also used the same procedure as 
iPorto de Mello et all (|2008l) . The Fe I lines span from 
to 5 eV in excitation potential, and from 10 to 115 mA 
in W\, allowing for an internally very precise solution of 
the excitation and ionization equilibrium. 

The formal standard error in the spectroscopic T e s was 
determined from the lcr uncertainty of the slope of the 
linear regression in the [Fe/H] vs. \ diagram, yield- 
ing the T c ff variation admissible at the lcr level. For 
the microturbulcnce velocity, the same procedure pro- 
vides the lcr microturbulence uncertainty in the [Fe/H] 
vs. W\ diagram. For the metallicity [Fe/H], we adopt 
the standard deviation of the distribution of abundances 
derived from the Fe I lines, which is larger than the 
errors in [Fe/H] due to T c g , £ and W\ errors. The 
error of the spectroscopic log g is estimated by evalu- 
ating the variation in this parameter which produces 
a disagreement of lcr between the abundances of Fe I 
and Fe II. The spectroscopic results thus determined for 
k 1 Cet are: T cff = 5780 ± 30 K, log 5 = 4.48 ± 0.10 dex, 
[Fe/H] = +0.07 ± 0.04 dex (for the Fe II lines, a = 0.05 
dex). The excitation & ionization equilibrium solution 
as a function of excitation potential is shown in Fig. [1] 
for all Fe I and Fe II lines. 

An additional T c g was determined by fitting the ob- 
served wings of Ha, followi ng closely the method of 
iLvra fc Porto de M ello (2005). This procedure is shown 
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Fig. 1. — The excitation & ionization equilibrium of Fe I and 
Fc II lines for re 1 Cet . 



in Fig . [2j Here we varied from ILvra fc Porto de Mellol 

(2005) routine by selecting manually individual flux 
points in the blue and red wings of the observed 
re 1 Cet Ha spectrum, free of telluric lines and photo- 
spheric line perturbations, in a total 21 and 43 points, 
respectively, for the blue and red wings. The T e g values 
derived from each wing agree within 5 K, and a mean 
value T e ff = 5645 ± 40 K. This error refers exclusively to 
the dispersion of T c s values attributed to the fitted profile 
data points. A m ore formal estimate of the unc ertainty 
of this procedure (|Lvra fc Porto de M ello 2005)), taking 
into account errors in the input atmospheric parameters 
and the continuum normalization, leads to an error of 50 
K. 

2.3. Effective temperature from the spectral energy 
distribution 

The photometric T e g- of k 1 Cet was calculated using a 
variety of methods. Firstly, we employed intermediate- 
band Stro mgren photometry coll ected from the GCPD 
database (Mermill iod et al.l 119971 ). which lists mean in- 
dices of V = 4.850 ± 0.008, (b - y) = 0.419 ± 0.003, 
mi = 0.235 ± 0.005, ci = 0.307 ± 0.003, and /3 = 
2.595 ± 0.016. We used the photometric grids of R. 
Napiwotzki (priv. comm.), which are based on ATLAS9 
model atmospheres. The determination of the T e ff was 
done by considering the lcr uncertainties in the color in- 
dices and following a Monte Carlo procedure. From 1000 
realizations we obtained a mean value and its correspond- 
ing standard deviation of T c s — 5650 ± 125 K. The er- 
ror bar only reflects the random uncertainty and does 
not account for systematic contributions. The Stromgren 
photometry also suggests roughly solar chemical compo- 
sition. 

A more accurate T e g determination was obt ained from 
the IRSED method of iMasana et all (|2006f l. This is 
based on the use of (2MASS) near-IR photometry and 
the fit of the spectral energy distributio n with stellar 
atmos phere models. The calibration of Masan a~et al.l 

(2006) takes into account possible systematic trends and 
applies corrections based on a sample of solar analogs. 
The photometry of k 1 Cet present in the 2MASS catalog 
is severely saturated but, fortunately, this star was se- 
lected as an ESQ standard and nu merous measurements 
are listed in iBouchet et al.l ()1991[ ). From the ESO IR 
photometry ( J ESO = 3.673±0.017, H ESO = 3.362±0.013, 
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Fig. 2. — Left: The blue wing of the Ha profile of k Cet, plotted with five representative spectra of the Moon, observed from 1995 to 
2002, and three theoretical models, centered at T e g = 5645 K and spaced by 50 K. The gray horizontal bars denote the positions of the 
stellar profile regions free from telluric and metal lines. Right: Same as above for the red wing. 



-Kf.ro = 3 . 282 ± 0.013), we applied transformations in 
ICarpenterl (|2001D to obtain photometry in the 2MASS 
system ( J 2M ass = 3.605±0.028, # 2 mass = 3.348±0.031, 
-ft^MASS = 3.239 ± 0.018). Using this photometry the 
IRSED method yields an effective temperature of T e g- = 
5685 ± 45 K when assuming solar values for logg and 
[Fe/H]. The dependence of the result on these adoptions 
is very weak, and k 1 Cet does have these parameters very 
similar to the Sun's, as already shown. 

2.4. Systematic offset between photometric, Ha and 
spectroscopic T e s values 

The three T e g values derived in the present work are 
not in agreement when the standard errors are consid- 
ered. The spectroscopic T c g is ~110 K higher than the 
Ha and photometric ones, which are in close agreement. 
When one considers the results of recent spectroscopic 
analyses of k 1 Cet, all of them based on high-quality, 
high- resolution spectra (Table [T]) , an interesting pattern 
emerges. Analyses based on photometric and B aimer 
line methods to derive T e s cluster very systematically 
in low values around T c ff = 5650 K, while those rely- 
ing on the excitation and ionization equilibria of Fe, or 
line depth ratios, cluster around T e g = 5750 K. A recent 
discussion on a possible offset between the spectroscopic 
T^ff scale on one side, and Balmer line and photometric 
ones on t he other, the former being th e hotter, has been 
given bv iPorto de Mello et all ([20081 ) . Authors do not 
generally agree whether this offset is due to NLTE ef- 
fects in cool stars, in the sense that for stars with T e g ~ 
5000 K the offsets are l arge, while for T eff ~ 6000 K 
good agreement is found (|Rarmrez et al.ll2007f) : this dis- 
crepancy is revealed in this case either by a disagree- 
ment between chemical abundances derived from differ- 
ent lines of the same species (atomic or molecular) or 
by the non- realization of the Fe I/Fe II ionization equi- 
librium. Another interpretation is the presence of high 
chromospheric activity, an enhanced non-local UV ra- 
diation field and a r esulting photospheric overionization 
(jSchuler et al.ll200"6T ). In a classical spectroscopic analy- 
sis, forcing agreement between Fe I and Fe II abundances 
under LTE, but in the presence of an overionizing radia- 
tion field, leads naturally to a higher T e g. As will be seen 
below, k 1 Cet is a very active star, with a chromospheric 
flux at the Ha core only slightly lower than a typical 



Hyades solar- type star, as well as a heightened UV flux 
with respect to the Sun below A1600 A, besides a very 
high X-ray luminosity. Since Fe II is essentially insen- 
sitive to NLTE effects for T e g va lues similar to that of 
the Sun (|Th evenin fe Idiartl[l999f) , the most reliable de- 
termination of the Fe abundance in this instance is that 
due to the Fe II lines, adopting the photometric and Ha 
profile T e ff. 

Accepting the likely presence of NLTE effects in k 1 Cet, 
we regard the following values as the most likely atmo- 
spheric parameters: T c g = 5665 ± 50 K, as a straight 
average of the T e g- values derived from photometry and 
the Ha, log 5 = 4.49 ± 0.10 and [Fe/H] = +0.10 ± 0.05 
dex, the latter due exclusively from the Fe II lines. Note 
that our metallicity value is in very good ageement with 
recent determinations (e.g.. iValenti fe Fischer|[2005l) . 

3. MAGNETIC ACTIVITY AND AGE 

k 1 Cet has always been recognized as a magnetically 
active star, with levels typical of a young solar analog. A 
space-based photom etric study using the MOST satellite 
( Wal keret al.ll2007h provided a precise value for its equa- 
torial rotation period of 8.77 days, which is some 3 times 
faster than that of today's Sun. This gives rise to a sig- 
nificantly enhanced magnetic dynamo and consequently 
stronger magnetically-generated phenomena. The mag- 
netic properties of n 1 Cet, together with spectroscopic 
diagnostics, can be used to place constraints on the age 
of the star. 

3.1. Ha absolute chromospheric flux 

Absolute chromosp heric radiative losses at t he co re of 
Ha were derived by iLvra fe Porto de Mellol (|2005h for 
a large sample of solar-type stars, including stars from 
Pleiades and Hyades clusters, and the Ursa Major kine- 
matic group. The average fluxes (and dispersions) of 
three Pleiades stars, seven UMa Group stars and seven 
Hyades stars are, respectively, 15.0 ± 1.9, 9.5 ± 2.9 and 
7.2 ± 1.2 (units are 10 5 erg cm" 2 s -1 A -1 ), in a clear 
age progression. ILvra fe Porto de Mellol (pOoBh estimate 
their Ha flux uncertainty (probably underestimated be- 
cause it does not account for errors in the photometric 
calibration of the fluxes) at 0.5 in the same units. In 
comparison to more classical spectroscopic indicators of 
chromospheric losses such as the H & K lines of Ca II 
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Fig. 3. — The Ha line core of re 1 Cet (thick full line), the Sun 
(thin full line), and mean spectra of representative stars of the 
Pleiades (dashed line) and Hyades (dotted line) clusters, and the 
Ursa Major moving group (dash-dotted line). 



(|Pasauinilll992l ). Ha fluxes are more subject to observa- 
tional uncertainty, but are also less sensitive to activity 
cycle phase and ro tational modulations, as well as tran- 
sient phenomena (|Lvra fc Porto de Mellol f2005f) . They 
should therefore be representative of the average level of 
stellar activity. 

In Fig. [31 mean spectra of the stellar groups, k 1 Cet and 
the Sun are overplotted. It is apparent that the mean 
chromospheric filling at the Ha core is much higher for 
the Pleiades, and similar for the UMa group and the 
Hyades. The appearance of the core flux spectra does 
not translate directly to the chromospheric fluxes, since 
the core flux is integrated in a 1.7 A- wide window and 
T e ff differences among the group members are not neg- 
ligible. Considering the dispersions of the mean Ha 
core fluxes of the Hyades cluster and the UMa group, 
nearly the same mean Ha activity level can be ascribed 
to them, and n 1 Cet is therefore compatible with their 
age range. Ag es for these stel l ar gro ups are 0.1 Gyr for 
the Plei ades dSchilbach et aTl|1995lh 0.63 Gyr for the 
Hyades dPerryman et al.l I1998T) and ~0.5 Gyr for the 
UMa Gro up (iKing et al.ll2003D. Applyi ng the flux cali- 
bration of lLvra fc Porto de Mellol (|2005l) to our Ha spec- 
trum leads to 7.3 • 10 5 erg cm~ 2 s _1 A -1 , which places 
k 1 Cet at the same flux level than an average Hyad, and 
below the flux level of an average member of the Ursa Ma- 
j or Gr oup. The age calibration of lLvra fc Porto de Mellol 
(2005) yields 0.65 Gyr for k 1 Cet, taking its chromo- 
spheric Ha flux at face value. 

3.2. Lithium abundance 

The lithium abundance of k 1 Cet was derived from the 
Li I resonance transition at A6707. A synthetic spec- 
trum was fitted to the FEROS spectrum, for two sets of 
atmospheric parameters: the purely spectroscopic solu- 
tion, T eff = 5780 K, log g = 4.48 dex, [Fe/H] = +0.07 
and £ = 1.21 km s , and the photometric/Ha solu- 
tion, T cff = 5665 K, log g = 4.49 dex, [Fe/H] = +0.10 
and £ = 1.20 km s . Mo del atmospheres we re inter- 
polated in the Kurucz grid (jKurucz et al.l[l993h and the 
synt hetic spectra w ere calculated with the MOOG rou- 
tine (|Sneden Ill973h . The FEROS instrumental broaden- 
ing profile was set at 0.07 A. The synthesis of Fe I lines 
in the A6707 vicinity provides a projected rotational ve- 



§ 0.9 
| 0.85 
3 0.8 



0.7 
0.08 
0.04 
' 
'-0.04 
-0.08 



(Lill 6707. a A) 



HD 20630 

Teff = 5780 logn(Li) = 2.f 2 



6702 6704 6706 6708 6710 6712 

Wavelength (A) 

Fig. 4. — Spectral synthesis of the A6707 Li I line of re 1 Cet: the 
purely spectroscopic set of atmospheric parameters was used. 

locity v sinz = 5± 1 km s _1 . For the spectroscopic atmo- 
spheric parameters, we derived log A?" (Li) = 2.12 (Fig. [4]), 
in the usual scale where log N (H) = 12.00. This value 
is in very goo d agre eme nt with the determin ations of 
iLuck fc Heiterl (|2006D and lPasquini et al.l (|1994l ). respec- 
tively, log AT (Li) = 2.04 and 2.13. 

This Li abundance places k 1 Cet, for the correspond- 
ing T G fj, in very good agreement with the Li sequence 
for the Praesepe and H yades clusters as determined by 
iSoderblom et a l. (1993). The two clusters are thought to 
be coeval at ~ 0.6 Gyr. Allowing for the spread in the 
Li abundances of the clusters' members, and the intrin- 
sically poor ability of Li abundances to discriminate age, 
this result does not actually constrain the age of k 1 Cet, 
but suggests that it is very young and not abnormal in 
its Li depletion history. 

3.3. Age 

Rotation periods for k 1 Cet have been estimated from 
spot modulation s. The determinations range between 
8.9 and 9.4 days (IGudel et al.|[l997l: iBaliunas et aT] |l995l: 
iMessina fc Guinanll2003l: iRucinski et al.ll2004D . The vari- 
ation is probably caused by the interplay of differential 
rotation and spots arising at different stellar latitudes 
of a star rotating w ith P = 8.77 days in the equator 
( Wal ker et al1l2007l) . In any case, the rotation period is 
within t he typical range for solar-type stars in the Hyades 
cluster (Radick et al. 1995), although close to the high 



end (IGudel et al.lll997l). The X-ray luminosity (log L x = 
28.8 in cgs: IGudel et al.lll997f) is also comparable to that 
of th e Hyades solar analogs (jBarrado v Navascues et al] 
1998). The same compatibility with Hyades members is 
also encountered in the case of chromospheric Ha emis- 
sion and in the abundance of Li in the atmosphere. All 
these indicators seem to suggest an age in the range 0.6- 
0.8 Gyr. In contrast, some age determinations based on 
chromospheric emission indices (like the Ca II H&K in- 
dex log seem to suggest a younger age around ~ 
0.4 G yr (|Lachaume et alj|1999t iMamaiek fc Hillenbrand! 
l2008h . 

All the age determination methods employed above 
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have significant uncertainties associated. For example, 
the rotation period vs. age relationship may carry a sig- 
nificant cosmic dispersion since not all stars have the 
same initial conditions and the same spin-down proper- 
ties. The Li abundance determination has been observed 
to have nearly one order of magnitude dispersion even 
within coeval clusters. And, finally, chronology with ac- 
tivity indicators may suffer also from cosmic dispersion, 
metallicity effects, etc. In summary, attributing a sin- 
gle age value for n 1 Cet is difficult and uncertain, but, 
putting all indicators together, we conclude that its age is 
likely to be in the range 0.4-0.8 Gyr (i.e., 0.6 ±0.2 Gyr), 
which is the value we adopt here. 

4. EVOLUTIONARY STATE 

Given the level of magnetic activity of k 1 Cet, it is ex- 
pected that the star still remains very close to, but not 
exactly on, its zero-age main sequence (ZAMS) position. 
We plotted k 1 Cet in two theoretical HR diag rams of 
the Yonsei-Yale (jKim et al.ll2002t lYi et alj|200a ) suite of 
models (Y 2 ). The first corresponds to the purely spec- 
troscopic solution of the atmospheric parameters (Fig. [5] 
top), Tee = 5780 K and [Fe/H] = +0.07, and the sec- 
ond to the photometric/Ha determination (Fig. [5] bot- 
tom), T e g = 5665 K, and the Fe abundance from Fe II 
lines, [Fe/H] = +0.10. The absolute magnitude from 
the Hipparcos parallax (109. 8 ± 0.78 mas) , coupled to 
the bolometric correction of iFlowerl (|1996l ) , results in 
logL/L = -0.080 ± 0.016, for T eB = 5780 K, and 
logi/L Q = -0.070 + 0.016, for T eB = 5665 K (both in a 
scale in which Mboi© = 4.75 and SCW© = —0.07). For 
the T e g , an uncertainty of 50 K was adopted. The dia- 
grams were slightly displaced in T e s and logi/L© (by 
+0.00186 and +0.0106 in log, respectively) so that a 
solar metallicity and solar mass track reaches the po- 
sition of the Sun exactly at T R ff = 5780 K for 4.53 Gyr 
(jGuenther fe Demarquelll997f) . This is acceptable for our 
current differential study. 

To study the evolutionary state of k 1 Cet in a more 
thorough manner we performed a series of Monte Carlo 
simulations taking into account both the "hot" and 
"cool" scenarios, with their corresponding luminosities 
and metallicities. For this purpose we have generated 
random values (N — 100000) with a Gaussian distri- 
bution for each of the input parameters, considering 
the errors as the standard deviations of the distribu- 
tions. For each combination of the values generated (T c ff, 
log(L/L©) and [Fe/H]) we carried out an interpolation 
in the synthetic Y 2 tracks to obtain the value of the mass 
and the age. Those parameter combinations that yielded 
unphysical situations (i.e., below the ZAMS) were dis- 
carded. In the case of the "cool" scenario, 85603 param- 
eter realizations could be used while only 19110 led to 
physical solutions in the case of the "hot" scenario. 

For our statistical study we only took into account 
those parameter combinations that resulted in an inter- 
polated age in the range 0.4-0.8 Gyr, which were 4047 
and 3715 for the "cool" and "hot" scenarios, respectively. 
According to this result there is no significantly bet- 
ter scenario. In both, we calculated the averages and 
standard deviations of the mass (the average of the age 
was obviously close to 0.6 Gyr), which were found to be 
1.045+0.011 M Q and 1.036+0.012 M Q for the "cool" and 
"hot" scenarios, respectively. In the case of the input pa- 



rameters, the average of the solutions yielding the correct 
age are T eff = 5705 ± 30 K, log(£/L Q ) = -0.076 + 0.014 
and [Fe/H] = +0.12 ± 0.05 for the "cool" scenario, and 
T off = 5740 ± 30 K, log(£/L ) = -0.073 ± 0.015 and 
[Fe/H] = +0.06 ± 0.05, for the "hot" scenario. The 
luminosities are essentially identical, while the tempera- 
tures are both pushed into mutual agreement, suggesting 
that a value close to the average is favored by the the- 
oretical models. Regarding the metallicity, the behav- 
ior is clearly different and the model calculations for the 
"cool" and "hot" scenarios tend to prefer slightly diver- 
gent values. This is expected since the metallicity has 
a strong effect in shifting the theoretical tracks, which, 
taking into account that the T e g and [Fe/H] uncertain- 
ties dominate the error budget, naturally constrains the 
physically acceptable solutions for each scenario within a 
narrow range of metallicity. Fig. [6] shows the histograms 
of the input and output parameters for the "cool" and 
"hot" scenarios. The solid lines depict the distribution 
from the entire simulation while dashed lines show the 
distribution of those solutions that yield an age in the 
interval 0.4-0.8 Gyr. 

The simulations favor slightly the "cool" spectroscopic 
solution and, although the statistical significance of this 
result is low, when considered with the likely presence 
of non-local radiative fields distorting the Fe I and Fe II 
populations, we find the "cool" spectroscopic solution 
more consistent and adopt it in the subsequent discus- 
sion. 

The position of k 1 Cet in the theoretical HR diagram, 
along with its magnetic activity context, is thus compat- 
ible with a slightly metal-rich star, slightly more massive 
than the Sun, and ~ 0.4-0.8 Gyr old. These T off and 
logL/L© values, along with M — 1.04 M©, when put 
into the well known equation: 

lead to log g — 4.49, in excellent agreement with the 
spectroscopic solution, and confirming that k 1 Cet is very 
close to the ZAMS. This surface gravity estimate is very 
insensitive to uncertainties in the inferred mass. 

Considering the Galactic orbit of k 1 Cet in the 
context of nearby s olar-t ype stars as analyzed by 
iPorto de Mello et all poM ). k 1 Cet is seen to have one 
of the lowest orbital eccentricities in the local population, 
and a mean galactocentric radius nearly identical to that 
of the Sun, and therefore a very similar Galactic orbit. 

5. UV IRRADIANCE 

Because of its apparent brightness, k 1 Cet has been 
subject to inten sive scrutiny with a variety of X-ray and 
UV telescopes (|Gudel et all \l99% iTelleschi et all |200H 
Paper I) that have revealed fluxes significantly higher 
than those of the current Sun. Paper I carried out an 
analysis of high-energy data with a wide wavelength cov- 
erage and found X-ray and EUV fluxes for k 1 Cet some 
20 and 10 times stronger, respectively, than today's Sun. 
In the FUV and UV, beyond the H Lyman a line (121.5 
nm) only fluxes for some strong f eatures have been pre- 
sented, mostly from IUE spectra ()Avreslll997l Paper I), 
but not the needed overall irradiance (strong lines + 
pseudo-continuum) . 
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Fig. 5. — Left: The state of evolution of re Cet in a suite of models with [Fe/H] = +0.07, the spectroscopic metallicity. The star is 
plotted with the spectroscopic T a g = 5780 K. Evolutionary tracks are labeled in solar masses. Note that the diagram does not correspond 
to solar metallicity. Right: Same as Fig[5]for [Fe/H] = +0.10. Dots along the tracks are labeled by ages in Gyr. Loci of same age between 
the 0.98 and 1.04 Mq tracks are linked by thin lines, re 1 Cet is plotted with the average T c ff from photometry and Ha profile fitting, T e s = 
5665 K. 



- 'Cool' scenario 

1 10000 wn 




Mass (M Q ) Age (Gyr) 



10000 r 




-0.10 0.00 
log(L/L ) 



5600 



5800 



T-„ (K) 



0.00 0.20 
[Fe/H] 
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0.4-0.8 Gyr. 
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We have compiled flux data of k 1 Cet in the UV and 
FUV. Spectra covering from 93 to 118 nm are available 
from FUSE and were discussed by iGuinan et al.l (2003) 
and in Paper I. Here we have used only night-time data 
since the relevant H Lyman lines are prone to contamina- 
tion from geocoronal emission. The interstellar medium- 
corre cted profile of the strong H Ly a line was taken 
from I Wood et al. (2005), which agrees well with the val- 
ues provided in Paper I. The rest of the spectrum (i.e., 
from H Lya up to nearly 300 nm) was taken from the 
CoolCAT catalog, which compiles echelle spectroscopy 
obtained with the HST STIS instrument. The combined 
spectrum, covering from 93 to 299 nm (with a few small 
gaps), was smoothed using a convolution with a gaussian 
of 0.3 nm FWHM and resampled at steps of 0.01 nm. The 
random uncertainties after such convolution are negligi- 
ble above 200 nm (well below 1%) and then slowly in- 
crease towards shorter wavelengths to reach about 10% 
at 140 nm and 20% at 93 nm. At long wavelengths the 
dominant error will be that of the standard flux calibra- 
tion but, according to STIS specifications and the infor- 
mation in CoolCAT, the systematic uncertainty should 
not exceed 5%. 

The k 1 Cet spectrum w as compared wi t h a sp ectrum 
of the Sun obtained from iThuillier et all (|2004f ). corre- 
sponding to medium solar activity. Fig. [7] illustrates such 
comparison. The left panel shows the normalized flux 
density at a distance of 1 AU from both k 1 Cet and the 
current Sun, while the right panel depicts the normalized 
flux ratio. The results clearly show that k 1 Cet is about 
35% fainter than today's Sun for wavelengths above 210 
nm, a range completely dominated by thermal radiation 
(except for the chromospheric emission of the Mg II h&k 
lines). Below 210 nm the difference in the ratio decreases 
to about 17% at 200 nm, 10% at 190 nm and the fluxes 
equal on average around 170 nm. Some chromospheric 
lines start to appear below 185 nm and those are much 
stronger in the case of k 1 Cet. Shortward of 170 nm the 
relative flux of k 1 Cet increases up to about a factor of 4 
at 120 nm and then even higher below. At those wave- 
lengths most of the flux is dominated by chromospheric 
lines that are significantly stronger for k 1 Cet. 

The results in Fig. [7]clearly illustrate how the Sun, be- 
cause of its higher temperature, possesses stronger photo- 
spheric UV radiation, yet, when the chromosphere emis- 
sion is considered, k 1 Cet emissions are stronger because 
of its higher level of magnetic activity. The transition 
between photosphere and chromosphere in terms of emis- 
sion is usually set at about 170 nm. However, note how 
in the interval between 170 nm and 210 nm, n 1 Cet is 
relatively brighter than expected. The flux ratio differ- 
ence with the Sun is about 10-15% while about 35% is 
fo und at longer wavelen gths . 

iCnossen et al.l (|2007f) presented estimates of the UV 
flux for k 1 Cet. The authors did not use measured UV 
data but employed astrophysical plasma models to gener- 
ate synthetic flux distributions from an emission measure 
distribution (EMD). Such EMD was obtained from the 
analysis of X-ray and EUV data, typically arising from 
coronal emission. The calculated UV fluxes, originating 
from much cooler plasma in the chromosphere, should 

8 http:/ /casa. colorado.edu/~ayres/CoolCAT/ 



thus be a crude appro ximation of rea l ity. C omparison of 
the fluxes obtained by ICnossen et al.l (|2007l ) is not possi- 
ble beyond a general apparent agreement, but real data, 
also with much higher resolution, ought to be preferred. 

6. A YOUNG SUN 

Evidence presented here for the young solar analog 
k 1 Cet suggests that the Sun was significantly more ac- 
tive in its past. This early activity for Sun- like stars 
is known in the X/EUV range and its evolution with 
the age of the star was studied in Paper I. Interestingly, 
k 1 Cet also exhibits enhanced fluxes in the ~ 100-200 
nm range. These emissions may have had an impact on 
the early evolution of Earth's atmosphere and may have 
played a role in the origin and development of life on 
Earth: it has been proposed, for instance, that UV radi- 
ation may help the synthesis of co mplex ribonucleotide s 
in plausible early Earth conditions (jPowner et alj |2009) . 
In this wavelength range, the solar emission drives the 
photochemistry and thus the molecular composition of 
planetary atmospheres. This can be gauged in Fig. HI 
where we show UV photoabsorption cross sections of rel- 
evant molecules. 

To illustrate how the spectral irradiance in the 100- 
300 nm domain can influence the photochemistry, we 
have computed the photolysis rates in an early Earth 
atmosphere subject to two different irradiance spectra. 
The first spectrum is that of k 1 Cet. The second spec- 
trum ( "model" ) is an educated guess of the UV spectral 
irradiance of the Sun at 0.6 Gyr based on our knowledge 
of stellar evolution and spectral synthesis, prior to the 
latests results of the Sun in Time program in Paper I, 
and on the basic assumption that the solar activity has 
not changed during the last 3.9 Gyr. This model of the 
young Sun spectrum is obtained by adding a synthetic 
photospheric irradiance spectrum and an activity com- 
ponent. The synthetic spectrum was computed by F. 
Castelli with the ATLAS9 model for the characteristics 
of the Sun at 0.6 Gyr (L = 0.73 L Q , T cff = 5680 K), 
as derived from the Y 2 stellar evolution models. The 
short wavelength emission associated with activity is as- 
sumed to be that of today's Sun. Such activity com- 
ponent can be calc ulated as the d i fferen ce between the 
observed spectrum IThuillier et al.l ((2004) and the pho- 
tospheric synthetic spectrunfl The difference, namely 
the activity component, becomes significant below about 
200 nm and dominates below 175 nm. Both spectral ir- 
radiances are scaled to 1 AU. The spectra of this simple 
young Sun model and k 1 Cet are compared in Fig. |H1 

We have considered two different atmospheric compo- 
sitions (see Fig. [101 top): 

• A CQ 2-rich atmosphere described in iKastind 
(1993). In this model the enhanced level of CO2 
provides enough greenhouse warming to compen- 
sate for the faint luminosity of the Sun 3.9 Gyr ago. 
The abundance of H2 is obtained for the current 
volcanic emission and assuming a diffusion-limited 
rate for the escape of hydrogen to space. This at- 
mosphere contains no organic species. 

• A more reduced atmosphere excerpted from 

9 ATLAS9 spectra of the present Sun are available at 
http: / /wwwuser. oat. ts.astro.it/castclli/sun.html 
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Fig. 7. — Left: Comparison of the observed UV spectra of k 1 Cet and the current Sun. Right: Ratio of the observed UV spectra of 
k Cet and the current Sun. 
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Fig. 8. — Photoabsorption cross sections of some molecules sus- 
pected to have been present in early Earth's atmosphere. 



iPavlov et all (|200l . that differs fr om the previous 
one by a significant level of CH4 and some pho- 
tochemically produced h ydrocarbons ( in pa rticu- 
lar C 2 H 2 and C 2 H 6 ). In iPavlov et ail (pOOl . the 
source of CH4 is assumed to be biogenic methano- 
genesis, but abiotic sources associate d with hy- 
drothermal activity are also plausibl e (jSleep et al.l 
[200l lAlbarede fc Blichert-Toftll2009l) . 

Both compositions include 1 bar of N2 as the main con- 
stituent. For both cases, our atmospheric thermal profile 
has been deduced from a present day terrestrial profile, 
assuming a constant stratospheric temperature equal to 
200 K between 14 and 60 km. This is of course not 
consistent with the detailed compositions but has little 
impact on t he photo l ysis r ates. Also, the detailed co mpo- 
sition from|Kastin3 |l993t ) and IPavlov et al.l (|2001l) were 
computed with a photochemical model that assumed dif- 
ferent UV fluxes. The purpose here is not to run a consis- 
tent photochemical model, which will be the done in fu- 
ture studies, but to illustrate how the photodissociation 
rates for these atmospheric compositions are sensitive to 
an enhanced UV irradiance. 

Photodissociation rates Ji{z) (s _1 ) at the altitude z of 
the different absorbing species i included in the model 




Young Sun model 



2000 

Wovelength (A) 



Fig. 9. — The observed ft 1 Cet spectrum vs. a simple model of 
the young Sun (see text). 

were computed in the range of wavelength [Ai, A2] as 



Hz) 




q itj (X)(x i (X)F(X,z)dX 



(1) 



which requires beforehand the implementation of their 
absorption cross-sections <Xi(A), of their different pho- 
todissociation pathways j each characterized by a quan- 
tum yield qtj (A) and of the incident stellar UV flux at 
every level in the atmosphere F(X, z) as well. This inci- 
dent stellar flux F(X, z) was calculated as a function of 
the diurnally averaged unattenuated stellar flux at the 
top of the atmosphere -Foe, (A) by considering only the 
integrated molecular absorption and Rayleigh scattering 
within molecular nitrogen N2. 



F(X,z) = F(X,oo)e~ T ^'+ diCf{X ' z) 



(2) 



The photolysis rates computed for both atmospheric 
compositions and both UV spectral irradiances are 
shown in Fig. [I0j These rates are found to be enhanced 
by a factor of 2-3 in the stratosphere and mesosphere 
for H2O and CO2 and by a factor of 4 for organic com- 
pounds. These photodissociations produce the radicals 
that activate the photochemical evolution, which is con- 
trolled by several hundreds of individual reactions, cou- 
pled with the vertical mixing and the condensation of 
some species. The resulting set of equations behaves 
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strongly non-linearly and can be very sensitive to the 
photolysis rates, which is why the UV spectrum has a 
critical influence. Therefore, considering k 1 Cet as a 
close analog of our Sun in its early times, such enhanced 
photodissociation rates were likely to trigger a peculiar 
atmospheric chemistry that we plan to investigate shortly 
by using a photochemical model of the primitive Earth's 
atmosphere. The higher H2O photolysis, associated with 
the intense heating of the thermosphere by the strong 
EUV irradiance, would also result in a higher escape rate 
of hydrogen to space, particularly important for the early 
evolution of Mars and Venus. The effect of the relatively 
strong UV fluxes is thus still significant for a 0.6 Gyr 
old Sun-like star, and is likely to play an even more im- 
portant role for younger stars, and for the atmospheric 
processes that occurred during the earliest stages of our 
planet. 

At wavelengths above 200 nm, the photospheric emis- 
sion dominates and the spectrum can be computed using 
a stellar atmosphere model. However, the inaccuracy of 
synthetic spectra can have a significant effect on the pho- 
todissociation rates. This inaccuracy has various origins 
as, for instance, the opacities, th e elemental abund ances, 
the limb darkening/brightening. Edvardsson (2008j) com- 
pared different models with observed stellar spectra, 
above 300 nm, and reported the occurrence of numer- 
ous patterns with widths of 0.5-1.5 nm, below 450 nm, 
with systematic uncertainties of about 10%. We noted 
even higher differences below 300 nm when comparing 
synthetic spectra obtained for the same stellar parame- 
ters but with different models. We can see in Fig. |H] that 
above 200 nm the synthetic model used for the young Sun 
fluctuates around that of k 1 Cet. This difference can be 
attributed to both the inaccuracy of the model and to 
the fact that k 1 Cet is a close but not exact replica of 
the Sun, having a slightly different mass and metallicity. 

Between 150 and 200 nm, the photon flux increases 
by more than two orders of magnitude. Which means 
that a 10% error on a 1 nm bin at 200 nm represent 10 
times more photons than all the flux in a 1 nm bin at 
150 nm. We can thus wonder about the influence of this 
error on the photodissociation rates. In the atmosphere 
models we considered, most of the species are photodisso- 
ciated at wavelengths below 200 nm and most absorption 
cross-sections decrease by orders of magnitude between 
150 and 200 nm. As a consequence, variations of the flux 
above 200 nm have a limited impact on the dissociation 
rates. For some species that absorb at A > 200 nm (like 
CO2 or hydrocarbons), photolysis rates can be slightly 
higher in the lowest atmosphere when computed with the 
young Sun model. But with the atmospheric composi- 
tions we considered, we checked that this effect remains 
negligible by comparing the rates from Fig. [10] with rates 
computed when both spectra are set to the same values 
at wavelengths above 200 nm. We should stress however 
that other atmospheric compositions, involving species 
that have strong absorption cross-sections at A > 200 nm 
(like for instance SO2 or O3), would be much more sen- 
sitive to this part of the spectrum. This demonstrates 
that using both the observed spectrum of a young Sun 
proxy and a theoretical model may be necessary for a full 
description of planetary atmospheres. 

7. SUMMARY AND CONCLUSIONS 



In this paper we have carried out an in-depth study of 
the bright, nearby solar analog k 1 Cet. Several methods 
have been used to estimate its effective temperature and 
chemical composition, yielding preferred values of T e g- = 
5665 ± 30 K, and [Fe/H] = +0.10±0.05. The systematic 
offset between T f{ values obtained from photometry /line 
profiles and the excitation/ionization Fe I and Fe II equi- 
libria is evidence of non-LTE effects, probably related to 
UV overionization due to the strong magnetic activity of 
this star. Adopting as the best atmospheric parameters 
the photometry and Ha profile T c g and the Fe II metal- 
licity, we have been able to set constraints to the stellar 
age, which should be between 0.4 and 0.8 Gyr. All the 
information gathered indicates that k 1 Cet is a star with 
nearly 1 solar mass in a relatively unevolved evolutionary 
stage. As such, it is an excellent match to the Sun as it 
was some 3.7-4.1 Gyr ago. 

The radiation from the young Sun must have played an 
essential role in shaping the atmospheres of the Solar Sys- 
tem planets. In particular, the UV flux is responsible for 
the photochemical processes in the atmosphere. We have 
been able to compile data, taken both with FUSE and 
HST, covering the entire UV and it shows that n 1 Cet's 
flux is some 35% lower than the current Sun's between 
210 and 300 nm, it matches the Sun's at 170 nm and 
increases to at least 2-7 times higher than the Sun's be- 
tween 110 and 140 nm. We have compared these fluxes 
with a "theoretical" young Sun estimated by adding the 
current chromospheric flux to a photospheric model with 
the correct radiative properties. We have used a photo- 
chemical model to calculate the photodissociation rates 
of the most relevant molecules in the assumed composi- 
tion of early Earth's atmosphere. The results indicate 
that such rates should have been several times higher 
than those resulting from a simplistic "theoretical" solar 
spectrum. 

Our calculations demonstrate that self-consistent plan- 
etary atmosphere calculations must account for the much 
stronger photodissociating radiation of the young Sun. 
The resulting chemistry could be significantly different 
from that commonly assumed. This is obviously very 
relevant at a significant point in the Solar System evo- 
lution, when life was gaining a secure foothold on Earth 
and Mars lost its liquid water inventory. 
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Fig. 10. — Top: The two ear ly Earth atm o spher es considered. Left: The "standard" atmosphere from iKastinel {1993). Right: A 
more reducing composition, from Pavlov ct al. (2001). Middle: The corresponding photodissociation rates of relevant molecular species 
computed for two different input solar spectra: re 1 Cet (solid lines) and a theoretical young Sun (dashed lines). Bottom: The percentage 
of enhancement in the photodissociation rates when adopting the spectral irradiance of re 1 Cet. 
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